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Abstract 

Skyrmion walls are topologically-nontrivial solutions of the Skyrme system which 
are periodic in two spatial directions. We report numerical investigations which show 
that solutions representing parallel multi- walls exist. The most stable configuration 
is that of the square AT- wall, which in the iV — » oo limit becomes the cubically- 
symmetric Skyrme crystal. There is also a solution resembling parallel hexagonal 
walls, but this is less stable. 
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The Skyrme system, originally introduced as a model of nucleons, is an archetypal 
(3+l)-dimensional classical field theory admitting topological soliton solutions. Much is 
known about various types of skyrmion solutions, for example: isolated skyrmions in R 3 , up 
to relatively high charge [H El E]; a triply-periodic 'Skyrme crystal' [U El El [2] ; a doubly- 
periodic 'Skyrme domain wall' (?]; and various types of singly-periodic 'Skyrme chains' 

The purpose of this note is to investigate static iV-wall solutions, ie. the N > 1 general- 
ization of the single- wall fields discussed in [7] . If one has two (or indeed N) well-separated 
parallel walls, then the force between them can be made attractive by a suitable relative 
orientation of the fields. So one expects there to be solutions representing N walls bound 
together, although a priori the walls might merge together to form a single wall. 

We investigate this by numerical minimization of the energy, and our main findings are 
as follows. There are two obvious shapes for a single wall, namely square and hexagonal, 
and it is known [7j that the latter has slightly lower energy than the former. If walls are 
allowed to attract, then they do not merge, but remain identifiable as separate parallel 
walls. There is a stable bound configuration representing two parallel hexagonal walls, but 
this is not the lowest-energy 2-wall state. For N > 2, the lowest-energy state consists of 
N parallel square walls (each one being a square array of half-skyrmions), and as N — > oo 
this approaches the skyrme crystal. 

The energy density of a static SU(2)-valued Skyrme field U{x^) on IR 3 is defined to be 

S := -~tr(L^) - -^tr([Lj, Lj] [L i: Lj]), (1) 

where Li = U~ 1 dU/dx t , and x^ = (x,y,z) are the spatial coordinates. In what follows, 
let us write U — $4 + i^jCj, where <jj are the Pauli matrices, and the real 4- vector field 
$ = ($ 1; $ 2 , $ 3 , $ 4 ) satisfies $ ■ $ = 1. 

In this note, we deal with configurations which resemble iV walls or sheets, each parallel 
to the xy-plane: so the field is periodic in x and y (with periods L x and L y respectively), 
and satisfies the boundary condition 

(2) 

For N = 1, and more generally for N odd, one has a domain wall which separates two 
vacuum regions, where $4 = 1 and $4 = — 1 respectively; for iV even, one has the same 
vacuum on both sides of the multi-layered sheet. In the asymptotic region \z\ ^> 1, the three 
fields $j are small; and they satisfy the Laplace equation, since the energy density reduces 
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to £ pa (<9j$j) 2 . Assuming (without loss of generality) that L y > L x , we see by separating 
variables that the leading behaviour as \z\ — > oo is typically ~ Csin (fiy) exp (— (J,\z\), 
where \i = 2ir/Ly. In particular, the fields approach their asymptotic values exponentially 
quickly, with a scale determined by the larger of L x and L y . 
The topological charge Q (over a single cell) is 



Q— Qdxdydz, (3) 

JT 2 xR 

where 

Q = 2^ £i ^ LiL ^ ^ 
is the topological charge density We claim that Q is an integer. If N is even, then (J2J) 
allows us to regard <E> as being defined, for topological purposes, on T 2 x S 1 ; and then Q 
equals the degree of If iV is odd, then it is not quite so obvious why Q is an integer, 
but it follows from the theorem in the Appendix of [8J. The energy E is defined to be 

1 f 

E := — — - / Sdxdydz, (5) 
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and it satisfies the usual Faddeev bound E > Q. 

In what follows, we describe iV-wall configurations which were found by numerical 
minimization of the energy functional E. We used a first-order finite-difference scheme 
for E, with the spatial points (x, y, z) being represented by a rectangular lattice having 
lattice spacing h, and we applied conjugate-gradient minimization. The lattice error in E 
goes like h 2 , and we extrapolated the finite-/?, results for both E and Q to h — 0. The 
extrapolated value of Q then gives a measure of the remaining error, which for the situations 
described below turns out to be less than 0.2%. The boundary condition (|2J) was modelled 
by imposing $ 4 = 1 at z = —L z /2 and $4 = (— 1)^ at z — L z /2. As remarked above, the 
walls are exponentially localized in z, and so as long as L z is taken to be large enough, 
there is no discernable finite-size effect; a value of L z = 10 + 2N turns out to be sufficient 
for this. In each case, we adjusted the periods L x and L y to their optimal size, meaning 
that the energy-per-cell is made as small as possible. Numerical minima were randomly 
perturbed and then re-minimized, as a test of their stability. As initial configurations 
we used the same sort of 'rational map ansatz' as in [7], involving a Weierstrass elliptic 
function of x + iy (the lemniscatic form to get square symmetry, and the equianharmonic 
form to get hexagonal symmetry), together with a suitable profile function f(z) satisfying 
f(-L z /2) = and f(L z /2) = Ntt. 
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N 




L 


1 


1.068 


4.25 


2 


1.053 


4.47 


3 


1.048 


4.54 


4 


1.046 


4.58 


5 


1.044 


4.61 



Table 1: Energy E and cell-size L for the square A/"- wall 



The results are consistent with the anticipated general principle that the lowest-energy 
configurations are arrays of half-skyrmions. For an AT- wall, we expect that each fundamen- 
tal cell will contain a multiple of 4A" half-skyrmions, and therefore its topological charge 
Q will be a multiple of 2 AT; this indeed turns out to be the case. As mentioned above, the 
walls do not merge, but retain their identity; the location of each wall can be determined 
by looking at the locus where $4 = 0. 

The simplest case to describe is the square one, with L y = L x = L; our results for 
1 < N < 5 are summarized in Table 1, which gives the energy-per-charge E, and the 
optimal value of L, for each N. Pictures of the N = 1 and N = 2 cases are presented in 
Figure 1, together with a plot of the energy data in Table 1. Let us first comment on the 
data. The normalized energy E of the square A^-wall is surprisingly close to having a 1/AT- 
dependence (although there is no obvious reason why this should be so); and extrapolating 
on this basis gives E ~ 1.039 in the A" — > 00 limit. This is very close to the energy of the 
(triply-periodic) skyrme crystal, a cubic array in which each fundamental cube contains 
eight half-skyrmions: its energy-per-charge, computed using the method described above, 
is E = 1.038. Further support for the claim that the square Af-wall tends to the skyrme 
crystal as A" — > 00 comes from looking at the symmetries of the field These include, 
for example, the translations 

X^X+\L X => ($i,$2,$3,$4) (-$l,-$2,$3,$4), 

y^y+±L y ($!, $ 2 , $3, $4) ^ ($i,-$ 2 ,-$3,$4), 

Z p ^Z p+1 => ($l,$2,$3,$4) >-> ($1, -$2,$3,-$4), 

where the third translation (in z) denotes moving from the p th wall to the (p + l) st wall. 
These are exactly the same as the translation symmetries of the skyrme crystal [2\. The 
values for the optimal cell-length L = L x = L y are consistent with their approaching 
L = 4.7 as N — > 00, this being the cell-size of the skyrme crystal (and similarly the 
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distance between each parallel pair of walls is approximately 4.7/2, as one would expect). 

Each 3-dimensional plot in Figure 1 is an isosurface of the energy density £, namely 
where S equals 0.6 times its maximum value. For the square case, the plots are over four 
fundamental cells. One clearly sees square arrays of half-skyrmions. Observe that, for 
N = 2, the half-skyrmions are aligned in the ^-direction; the same is true for N > 2. 

Let us turn now to the case of hexagonal symmetry. For ease of computation, we 
follow the same scheme as in [7j, namely taking L y = \/3L x and fitting two fundamental 
parallelograms into the corresponding rectangle. Each such rectangle, of each wall, contains 
eight half-skyrmions, as is seen in the hexagonal 2- wall picture of Figure 1. For N = 1, 
the energy of the hexagonal arrangement is E « 1.062, less than that of the corresponding 
square case [7j; but for iV > 2, the hexagonal arrangement is less efficient than the square 
one, and (depending on the values of L x and L y ) it is either a local minimum of the energy 
functional, or it is unstable. There is a local-minimum hexagonal 2-wall solution with 
energy E « 1.055, which is only very slightly (less than 0.2%) higher than that of the 
square 2-wall. Its energy density is depicted in Figure 1; one feature to note is that the two 
walls are not aligned in the z-direction, but are offset. If L x and L y are allowed to change 
so that the relation L y = y3L x no longer holds, then this solution becomes unstable, and 
changes into the square 2-wall. 

An isolated skyrmion of charge Q > 3 typically has a polyhedral-shell structure, anal- 
ogous to carbon fullerenes, and it may be viewed as constructed from a section of the 
hexagonal 1-wall (graphene), with the insertion of defects to create a spherical shell [HE]. 
There has also been an investigation [9] of the possibility of constructing skyrmions as 
multi- walled spherical shells, with the 'shell material' consisting of a double or triple wall. 
For the cases that were examined in [9J, either the walls coalesced, or one obtained a 
structure which resembled a shell-like part of the skyrme crystal. The findings reported 
above are consistent with this; in particular, multiple hexagonal walls appear to be rather 
unstable, and therefore unsuitable for constructing shells. But it does raise the possibil- 
ity of stable high-charge skyrmions constructed as shells of square multi-wall material, or 
equivalently as hollow chunks of skyrme crystal, and this would be worth investigating 
further. 
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Figure 1: Energy densities of the square 1-wall, square 2- wall, and hexagonal 2-wall; and 
plot of the energy E for the square iV-wall (1 < iV < 5) and hexagonal iV-wall (1 < N < 2). 
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